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Introduction 


This report summarizes work accomplished under NASA Grant NCC2-5068, June 1 5, 
1994 to September 30, 1995, entitled "Advanced Methods of Structural and Trajectory Analysis 
for Transport Aircraft," The effort was in two areas (1) development of advanced methods of 
structural weight estimation, and (2) development of advanced methods of trajectory 
optimization. The majority of the effort was spent in the structural weight area. During the 
course of the grant, there were slight deviations from the original work statement due to changing 
priorities of the sponsor. M. Chambers and H.-C. Chou were the graduate student research 
assistants assigned to the project. 

Structural Weight Estimation 

Analytically-based weight estimating subroutines for the fuselage and wing structures of 
transport aircraft were developed and integrated into the ACSYNT vehicle synthesis code. The 
subroutines were based on previously developed codes for hypersonic aircraft [Refs 1-3]. 

The methods of analysis used to develop the weight estimating routines are discussed in 
detail in a NASA TM, attached as an Appendix to this report. The TM also serves as a users 
manual for running the routines as an integral part of ACSYNT. The methods used and the 
resulting program will be only briefly reviewed here. 

Preliminary weight estimates of aircraft traditionally have been made using empirical 
methods based on the actual weights of existing aircraft. The fuselage and wing designs of 
advanced aircraft, however, may be significantly different from those of existing aircraft. This 
means that an empirically based method of weight estimation may not be valid. On the other 
hand, finite-element methods of structural analysis, commonly used in aircraft detailed design, are 
not appropriate for conceptual and preliminary design, because of the large number of specific 
cases that need to be considered, and the large amount of input necessary for each case. 

The body structural weight estimation method developed for ASCYNT is based on a third 
approach, beam theory structural analysis. This results in a weight estimate that is directly driven 
by material properties, load conditions, and vehicle size and shape, and is not confined to an 
existing data base. Since the analysis is done station-by-station along the vehicle longitudinal axis 
and along the wing structural chords, the distribution of loads and vehicle geometry is accounted 
for, giving an integrated weight that accounts for local conditions. Because of the beam 
assumption, the analysis is valid only for high aspect ratio wings and high fineness ratio fuselages. 

Although the weight estimating routines are based on previously developed codes, there 
was substantial modification to the codes. Along the more significant changes were: (1) 
modeling the fuselage shell as a cylinder with two power-low ends (the previous program had two 
power low sections back-to-back); (2) changing the wing lift distribution from trapezoidal to 
elliptical; (3) allowing for linearly varying wing thickness ratio along the span; and (4) allowing 
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the engines to be mounted on either the fuselage or the wing, or a combination of both. The 
routines were integrated with ACSYNT and verified with test cases. The code was made more 
user friendly by streamlining and adding comment statements. 

The structural weight routines were used to estimate the fuselage and wing weights of 
nine existing transport aircraft. This was necessary to assess the validity of the methods and to 
determine the "nonoptimum" weight. The nonoptimum weight is the portion of the weight not 
estimated by the structural analysis, such as fasteners, doublers, cutout reinforcement, uniform 
gage penalties, etc. The results show that the developed weight routines give statistically very 
good weight estimates. 

Trajectory Optimization 

In the trajectory area, the main task was to add a trajectory optimization routine to 
ACSYNT. The optimization is based on the energy-state dynamic model, and can be used to 
minimize time, fuel, or a weighted combination of the two. 

The energy-state approximation has been used successfully many times to obtain guidance 
laws for a wide variety of aircraft and missions [Refs. 4-8], In particular, in Ref. 5 this method is 
used to derive optimal guidance laws for minimizing direct operating costs (DOC) for transport 
aircraft. The guidance laws were found to be extremely accurate and have been implemented for 
on-board control of existing airplanes 

The starting point of the analysis is the equations of motion of an aircraft center of mass 
flying in a vertical plane above a flat, non-rotating earth: 


h = siny 

v - 7c° s Y -D sin y 

m 

7sin y+Z - mg cos y 
mv 

m -- -CT 


( 1 ) 


where the state variables are altitude, ( h ), speed (v), flight path angle (y), and mass (/??); the forces 
are thrust (7), lift ( L ), and drag (. D ); a, g and C are the angle of attack, the gravitational 
acceleration, and the specific fuel consumption, respectively. Next introduce a new state variable, 
the total mechanical energy per unit weight, 

£ h* — v 2 (2) 

2 g 
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Taking the time derivative of this, using eqns. (1 ), and replacing v by E as a state variable gives 


E . —(Tcoea-D) - P 
f»g 

h - vsiny 

7sin a*L - mg cos y 
mv 

m - -CT 


( 3 ) 


where P is called the specific excess power. 

The energy-state approximation consists of assuming that both y and y are small, and thus 
that h is small as well (these are very good approximations for transport aircraft flight). 
Neglecting these terms gives 

E - P 

m - CT (4) 

o = Tsma.*L -mg 

which are the equations of the energy-state approximation. 

The quantity to be minimized along the trajectory is a weighted sum of flight time and fuel 
mass consumed: 


4> - K x uKjn f (5) 

This allows minimum time trajectories to be computed = 1, K 7 = 0 ), minimum fuel 
consumption (^T, = 0, K 2 = 1 ), and minimum DOC (by proper weighting of K x and K 2 ). For a 
given energy gain along the trajectory, the quantity to be minimized is 

♦/ V E f : 

J x = fd< f> = ji ult - fZ-dE (6) 

. *0 'o E o 


where eqn. (4) was used. It is assumed that <j> > , P > 0 and that E is monotonically increasing. 
For convenience, the integrand in eqn. (6) is inverted and we maximize 
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( 7 ) 


J - f—dE - (FdE 

i* { 

From eqns. (4), (5), and (7), the function F is: 

F v(7cos a-D) 

" mg(K^KfT) 

The optimal trajectory guidance algorithm is then 

v = arg max [F] E (9) 

v 

that is, at each energy level along the flight path, choose the v that maximizes F, subject to any 
relevant constraints. Once this v has been determined, the integration to the next step is done 
with eqns. (1) [or equivalently eqns. (3)], with two choices of approximation. One choice is with 
Y = 0 , and the other is with both y = 0 and h = 0 . The latter choice, corresponding to eqns. 

(4), is available in case severe jumps are present in the energy-state optimal path, making 
integration including the h term impossible. (The jumps which may occur are an important topic 
for future investigation.) 

The first example of use of the trajectory optimization was the Boeing 747. Figure 1 
shows the minimum time climb trajectory for this aircraft. Also shown are lines of constant E and 
the contours of the function F, illustrating that F is maximized at each value of E. Figure 2 shows 
the minimum fuel trajectory, and Figure 3 compares the two flight paths. It is seen that the 
minimum fuel path is somewhat higher than the minimum time path, the minimum DOC path 
would be between these two. Note that these paths have an acceleration at ground level to about 
M 0.5, which would not be operationally acceptable. All integrations were done with the h term 
included. 

The algorithm was also used to determine optimal climb trajectories of two high speed 
(supersonic) civil transport designs, designated here as HSCT1 and HSCT 2. Figure 4 shows the 
minimum time path for HSCT1 . Also shown are the contours of F (labeled P,/B). The maximum 
dynamic pressure constraint (q max) has been set unrealistically high in order to better illustrate 
the nature of the path. The trajectory follows the terrain and q max limits, except for a transonic 
dive and a brief low supersonic jump to lower q. Figure 5 magnifies this behavior in the 
transonic/low supersonic region. The "trough" in the performance function F(M, h), which causes 
this behavior, is clearly shown in Figure 6. 


4 


Figures 7, 8, and 9 show the same information for HSCT 1 for minimum fuel. As 
compared with minimum time paths, the minimum fuel paths have more jumps and are at lower q. 
Figure 9 shows that F(M,h) is very "hilly" in the transonic region for this case, having multiple 
local maximums. This makes searching for the global maximum very difficult. 

The elapsed time and weight of HSCT 1 as a function of E are shown in Figures 10 and 11. 
It is seen that the minimum time path takes less time and has more fuel expenditure than the 
minimum fuel path, thus validating the calculations As before, the minimum DOC path would lie 
between the two and would give intermediate time and weight. 

The minimum time and minimum fuel trajectories for HSCT2 are displayed in Figures 12 - 
15. The paths tend to follow the q max constraint, except for a transonic climb and dive; the dive 
is more severe for the minimum fuel path. 

In order to execute the trajectory optimization routines in HAVOC, the following 
variables need to be set: 


VIND 

set to -2 for climb trajectory optimization 

HDOTOP 

set to TRUE.:/j included in integration 


FALSE.: h not included in integration 

NSTEP 

number of search steps at each energy level (currently set to 
NELGCL, number of energy levels) 

WGTFACK1 

weighting parameter on fuel weight (also called FMNFUL) 

WGTFACK2 

weighting parameter on time (also called FMNTIM) 
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Figure 1 . The optimized flight path of B747 in Min. Time 
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Figure 2. The optimized flight path of B747 in Min. Fuel. 
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Figure 3. The optimized flight paths of B747. 
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Figure 4. The optimized flight path of HSCT1 in Min. Time. 
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Figure 5. The optimized flight path, in min Time, and cost functional contours of 

HSCT 1 in the transonic region. 
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Figure 6. Cost functional mesh and ascent trajectory of HSCT1 in Min. Time. 
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Figure 7. The optimized flight path of HSCT1 in Min. Fuel. 
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Figure 8. The optimized flight path, in Min. Fuel, and cost functional contours of 

HSCT1 in the transonic region. 
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Figure 10. Time consumptions vs. energy level of HSCT1. 
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Figure 11. Weight status vs. energy level of HSCT1. 



Alt. (ft) 



Mach No. 


Figure 12. The optimized flight path of HSCT2 in Min. Time. 
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Figure 13. Cost functional mesh and ascent trajectory of HSCT2 in Min. Time. 
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Figure 14. The optimized flight path of HSCT2 in Min. Fuel. 
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SUMMARY 


A method of estimating the load-bearing fuselage weight and wing weight of 
transport aircraft based on fundamental structural principals has been developed. This 
method of weight estimation represents a compromise between the rapid assessment of 
component weight using empirical methods based on actual weights of existing aircraft, 
and detailed, but time-consuming, analysis using the finite element method. The method 
was applied to nine existing subsonic transports for validation and correlation. Integration 
of the resulting computer program, PDCYL, has been made into the weights-calculating 
module of the ACSYNT computer program. ACSYNT has traditionally used only 
empirical weight estimation methods; PDCYL adds to ACSYNT a rapid, accurate means 
of assessing the fuselage and wing weights of unconventional aircraft. PDCYL also 
allows flexibility in the choice of structural concept, as well as a direct means of 
determining the impact of advanced materials on structural weight. 

Using statistical analysis techniques, a relation between the load-bearing fuselage 
and wing weights calculated by PDCYL and corresponding actual weights was 
determined. A User’s Manual, as well two sample outputs, one for a typical transport and 
another for an advanced concept vehicle, are given in the Appendices. 


A fuselage cross-sectional area 

A b fuselage surface area 

A F frame cross-sectional area 

( AR ) aspect ratio of wing 

b wing span; intercept of regression line 

b s stiffener spacing 

b s structural semispan of wing, measured along quarter chord from fuselage 
b w stiffener depth 

C F Shanley’s constant 

C P center of pressure 

C R theoretical root chord of wing 

C R root chord of wing at fuselage intersection 

C 5| portion of wing leading edge not used for structural box 


C $2 portion of wing trailing edge not used for structural box 

Csr structural root chord of wing 

Or structural tip chord of wing 

C T tip chord length of wing 

d frame spacing 

d w optimum web spacing of wing 

D maximum diameter of fuselage 

e wing buckling exponent 

e c wing cover material factor 

E Young’s modulus of shell material 

E F Young’s modulus of frame material 

Ecy compressive yield strength 

Fm ultimate tensile strength 

F s shear stress 

h thickness of sandwich shell 

h e step function for 1 th engine on wing 


h lgi step function for I th landing gear on wing 
I F frame cross-section moment of inertia 

I y moment of inertia about y-axis 

Iy 1ft s 

K fi frame stiffness coefficient, IjJ A F 

K mg shell minimum gage factor 

K p shell geometry factor for hoop stress 
K s constant for shear stress in wing 

K th sandwich thickness parameter 

l B fuselage length 

Ile length from leading edge to structural box at theoretical root chord 

1 TE length from trailing edge to structural box at theoretical root chord 

l NG length from nose to nose gear 

l MC length from nose to fuselage mounted main gear 

/ , nose length of fuselage 

l 2 tail length of fuselage 

/„ distance to breakpoint of fuselage 

L lift 

L T maximum vertical tail lift 

M longitudinal bending moment; wing longitudinal bending moment 
m buckling equation exponent; slope of regression line 
n normal load factor 

N z axial stress resultant 

N Xb bending stress resultant 

N Xp pressure stress resultant 

fJ x tensile axial stress resultant 

N~ compressive axial stress resultant 

N y hoop direction stress resultant 

n x longitudinal acceleration 

P perimeter 


P g internal gage pressure 

P s perimeter of shell 

P w perimeter of walls 

Pj exponent of power law of nose of fuselage 

P 2 exponent of power law of tail of fuselage 

r radius of fuselage 

r(y) total wing chord as a function of position along quarter chord 

r/y j structural wing chord as a function of position along quarter chord 

R correlation coefficient used for regression 

Pfm fineness ratio 

Rffj ratio of horizontal tail station to fuselage length 

R le ratio of leading edge station of wing to fuselage length 

R P ratio of length to leading edge of fuselage propulsion to fuselage length 

Rp 2 ratio of length to trailing edge of fuselage propulsion to fuselage length 

R ng ratio of length to nose gear to fuselage length 

R mg ratio of length to fuselage mounted main gear to fuselage length 

R,(y) thickness ratio of wing, linearly interpolated between root and tip thickness ratios 

R tap taper ratio of wing 

S B plan area of the fuselage 

Slg stroke of landing gear 

S P plan area of wing 

t(y) thickness of wing box as a function of position along quarter chord 

t c core thickness 

tf face sheet thickness 

t s skin thickness 

t g material gage thickness, t s /K mg 

t w stiffener thickness 

t„ g material minimum gage thickness 

t total equivalent isotropic thickness of shell and frames 

t B total equivalent isotropic thickness of fuselage structure 

t F smeared equivalent isotropic thickness of frames 


t s equivalent isotropic thickness of shell 

i s shell thickness required to preclude buckling failure 

t Sc shell thickness to preclude compressive failure 

t Sc shell thickness required to meet minimum gage constraint 

i Sj . shell thickness required to preclude tensile failure 

t T smeared tension tie thickness 

/ M , smeared wall thickness 

/ thickness of wall to meet minimum gage constraint 
t Wr thickness of wall required to prevent tensile failure 

T torque on wing carrythrough structure 

V B fuselage volume 

V w volume of wing structural box, including structure 
V ] volume of fuselage nose 

V 2 volume of fuselage tail 

w c width of carrythrough structure of wing (usually same as D) 

W body structural weight 

W weight of wing per unit span 

W, ideal fuselage structural weight 

Wfj fuel weight 

W NO weight of nonoptimum material 

W s vehicle longitudinal weight distribution 

W T0 gross takeoff weight 

W/S shell structural weight per unit surface area 

x longitudinal fuselage coordinate 

x eaic weight calculated by PDCYL 

x m distance from nose to theoretical quarter chord of horizontal tail 

x LE distance from nose to leading edge of wing 

x Pi distance from nose to leading edge of fuselage mounted propulsion 

x P) distance from nose to trailing edge of fuselage mounted propulsion 

y transverse fuselage coordinate; wing coordinate measured along 1/4 chord 


y ac , actual weight 

y est estimated weight after regression 

2 vertical body coordinate 

Z(y) total width of wing box as a function of position along quarter chord 
%s(y) width of wing box structure as a function of position along quarter chord 
5 frame deflection 

e shell buckling efficiency 

e c wing cover structural efficiency 

z w wing web structural efficiency 

A wing sweep 

p shell structural material density; wing structural material density 

p g gross fuselage density 

p F frame structural material density 

p wing loading 

ct 5 allowable shear stress for wing 

I sum over fuselage or wing length; solidity of wing 

\|/ truss core angle 


INTRODUCTION 


A methodology based on fundamental structural principles has been developed to 
estimate the load-carrying weight of the fuselage and basic box weight of the wing for 
aircraft, and has been incorporated into the AirCraft SYNThesis Program (ACSYNT). 
This weight routine is also available to run independently of ACSYNT, and is a 
modification of a collection of previously developed structural programs (References 1- 
4). The main subroutine called by ACSYNT is PDCYL. This study has concentrated on 
modem transport aircraft because of the detailed weight information available, allowing 
the weights output from PDCYL to be compared to actual structural weights. The 
detailed weight statements also allow non-optimum factors to be computed which, when 
multiplied by the load-bearing structural weights calculated by PDCYL, will give good 
representative total structure weight estimates. These non-optimum factors will be 
computed through a regression analysis of a group of nine transport aircraft. 

PDCYL is able to model both skin-stringer-ffame and composite sandwich shell 
fuselage and wing box constructions. Numerous modifications were made to PDCYL 
and its associated collection of subroutines. These modifications include the addition of 
detailed fuselage shell geometry calculations; optional integration of a cylindrical 
fuselage midsection between the nose and tail sections; addition of landing and bump 
maneuvers to the load cases sizing the fuselage; ability to introduce an elliptical 
spanwise lift load distribution on the wing; variation of wing thickness ratio from tip to 
root; ability to place landing gear on the wing to relieve spanwise bending loads; 
distribution of propulsion system components between wing and fuselage; and the 
determination of maximum wingtip deflection. 

Brief description of ACSYNT 

The Aircraft Synthesis Computer program, ACSYNT, is an integrated design tool 
used in the modeling of advanced aircraft for conceptual design studies. ACSYNT 
development began at NASA- Ames Research Center in the 1970's and continues to this 
day. The ACSYNT program is quite flexible and can model a wide range of aircraft 
configurations and sizes, from remotely piloted high altitude craft to the largest transport. 

The ACSYNT program uses the following modules, not necessarily in this order: 
Geometry, Trajectory, Aerodynamics, Propulsion, Stability, Weights, Cost, Advanced 
Aerodynamic Methods, and Take-off. An ACSYNT run would normally progress as 
follows: the Geometry module is called to define the aircraft shape and configuration; 
the Trajectory module then runs the vehicle through a specified mission; finally the 
Weight and Cost modules are executed. To determine the performance of the vehicle at 
each mission point, the Trajectory module will call the Aerodynamics and Propulsion 
modules. 

After the mission is completed, the calculated weight of the aircraft may be 
compared with the initial estimate and an iteration scheme run to converge upon the 


appropriate aircraft weight. This process is necessarily iterative as the aircraft weight 
ACSYNT calculates is dependent upon the initial weight estimate. 

ACSYNT is able to perform a sensitivity analysis on any design variable, such as 
aspect ratio, thickness-to-chord ratio, fuselage length or maximum fuselage diameter. 
Sensitivity is defined as (change in objective function/value of objective function) 
divided by (change in design variable/design variable). As an example, if gross weight is 
the objective function and decreases when the wing thickness-to-chord ratio increases, 
then the sensitivity of thickness-to-chord ratio is negative. It is important to note that 
while this increase in thickness-to-chord ratio lowers the gross weight of the aircraft, it 
may also have a detrimental effect on aircraft performance, possibly violating 
performance constraints. 

ACSYNT is also able to size multiple design variables by optimizing the objective 
function, subject to some preset limits or constraints. Objective functions such as gross 
weight, performance, or fuel weight, may be maximized or minimized subject to some 
upper or lower bounds. During the optimization process, lower and upper bounds may 
also be imposed on the design variables. 

Three Methods of Weight Estimation 

Two methods are commonly available to estimate both the load-bearing fuselage 
weight and wing box structure weight of aircraft. These methods, in increasing order of 
complexity and accuracy, are empirical regression and detailed finite element structural 
analysis. Each method has particular advantages and limitations which will be briefly 
discussed in the following sections. There is an additional method based on classical 
plate theory (CPT) which may be used to estimate the weight of the wing box structure. 

Empirical 

The empirical approach is the simplest weight estimation tool. It requires 
knowledge of fuselage and wing weights from a number of similar existing aircraft in 
addition to various key configuration parameters of these aircraft in order to produce a 
linear regression. This regression is a function of the configuration parameters of the 
existing aircraft and is then scaled to give an estimate of fuselage and wing weights for an 
aircraft under investigation. Obviously, the accuracy of this method is dependent upon 
the quality and quantity of data available for existing aircraft. Also, the accuracy of the 
estimation will depend on how closely the existing aircraft match the configuration and 
weight of the aircraft under investigation. All of the empirical regression functions 
currently in the ACSYNT program give total fuselage weight and total wing weight. 

Finite Element 

Finite element analysis is the matrix method of solution of a discretized model of 
a structure. This structure, such as an aircraft fuselage or wing, is modeled as a system of 
elements connected to adjacent elements at nodal points. An element is a discrete (or 
finite) structure that has a certain geometric makeup and set of physical characteristics. A 
nodal force acts at each nodal point, which is capable of displacement. A set of 


mathematical equations may be written for each element relating its nodal displacements 
to the corresponding nodal forces. For skeletal structures, such as those composed of 
rods or beams, the determination of element sizing and corresponding nodal positioning 
is relatively straightforward. Placement of nodal points on these simple structures would 
naturally fall on positions of concentrated external force application or joints, where 
discontinuities in local displacement occur. 

Continuum structures, such as an aircraft fuselage or wing, which would use some 
combination of solid, flat plate, or shell elements, are not as easily discretizable. An 
approximate mesh of elements must be made to model these structures. In effect, an 
idealized model of the structure is made , where the element selection and sizing is 
tailored to local loading and stress conditions. 

The assembly of elements representing the entire structure is a large set of 
simultaneous equations that, when combined with the loading condition and physical 
constraints, can be solved to find the unknown nodal forces and displacements. The 
nodal forces and displacements are then substituted back into the each element to produce 
stress and strain distributions for the entire structural model. 

Classical Plate Theory' 

CPT has been applied to wing structure design and weight estimation for the past 
20 years. Using CPT a mathematical model of the wing based on an equivalent plate 
representation is combined with global Ritz analysis techniques to study the structural 
response of the wing. An equivalent plate model does not require detailed structural 
design data as required for finite element analysis model generation and has been shown 
to be a reliable model for low aspect ratio fighter wings. Generally, CPT will 
overestimate the stiffness of more flexible, higher aspect ratio wings, such as those 
employed on modem transport aircraft. Recently, transverse shear deformation has been 
included in equivalent plate models to account for this added flexibility. This new 
technique has been shown to give good representations of tip deflection and natural 
frequencies of higher aspect ratio wings. No fuselage weight estimation technique which 
corresponds to the equivalent plate model for wing structures is available. 

Need for better, intermediate method 

Preliminary weight estimates of aircraft are traditionally made using empirical 
methods based on the weights of existing aircraft, as has been described. These methods, 
however, are undesirable for studies of unconventional aircraft concepts for two reasons. 
Hirst, since the weight estimating formulas are based on existing aircraft, their application 
to unconventional configurations (i.e., canard aircraft or area ruled bodies) is suspect. 
Secondly, they provide no straightforward method to assess the impact of advanced 
technologies and materials (i.e., bonded construction and advanced composite laminates). 

On the other hand, finite-element based methods of structural analysis, commonly 
used in aircraft detailed design, are not appropriate for conceptual and preliminary design, 
as the idealized structural model must be built off-line. The solution of even a 


moderately complex model is also computationally intensive and will become a 
bottleneck in the vehicle synthesis. Two approaches which may simplify finite-element 
structural analysis also have drawbacks. The first is to create detailed analyses at a few 
critical locations on the fuselage and wing, then extrapolate the results to the entire 
aircraft, but this can be misleading because of the great variety of structural, load, and 
geometric characteristics in a typical design. The second method is to create an 
extremely coarse model of the aircraft, but this scheme may miss key loading and stress 
concentrations in addition to suffering from the problems associated with a number of 
detailed analyses. 

The fuselage and wing structural weight estimation method employed in PDCYL 
is based on a third approach, beam theory structural analysis. This results in a weight 
estimate that is directly driven by material properties, load conditions, and vehicle size 
and shape, and is not confined to an existing data base. Since the analysis is done station- 
by-station along the vehicle longitudinal axis, and along the wing structural chord, the 
distribution of loads and vehicle geometry is accounted for, giving an integrated weight 
that accounts for local conditions. An analysis based solely on fundamental principles 
will give an accurate estimate of structural weight only. Weights for fuselage and wing 
secondary structure, including control surfaces and leading and trailing edges, and some 
items from the primary structure, such as doublers, cutouts, and fasteners, must be 
estimated from examination of existing aircraft. 

The equivalent plate representation, which is unable to model the fuselage 
structure, is not used in PDCYL. 


METHODS 


Overview 

Since it is necessary in systems analysis studies to be able to rapidly evaluate a 
large number of specific designs, the methods employed in PDCYL are based on idealized 
vehicle models and simplified structural analysis. The analyses of the fuselage and wing 
structures are performed in different routines within PDCYL, and, as such, will be 
discussed separately. The PDCYL weight analysis program is initiated at the point where 
ACSYNT performs its fuselage weight calculation. PDCYL first performs a basic 
geometrical sizing of the aircraft in which the overall dimensions of the aircraft are 
determined and the propulsion system, landing gear, wing and lifting surfaces are placed. 

Fuselage 

The detailed fuselage analysis starts with a calculation of vehicle loads on a 
station-by-station basis; Three types of loads are considered - longitudinal acceleration 
(applicable to high-thrust propulsion systems), tank or internal cabin pressure, and 
longitudinal bending moment. All of these loads occur simultaneously representing a 
critical loading condition. For longitudinal acceleration, longitudinal stress resultants 
caused by acceleration are computed as a function of longitudinal fuselage station, these 
stress resultants are compressive ahead of the propulsion system and tensile behind the 
propulsion system. For pressure loads, the longitudinal distribution of longitudinal and 
circumferential (hoop) stress resultants is computed for a given shell gage pressure 
(generally 12 psig). There is an option to either use the pressure loads to reduce the 
compressive loads from other sources or not to do this, in either case, the pressure loads 
are added to the other tensile loads. 

Longitudinal bending moment distributions from three load cases are examined for 
the fuselage. Loads on the fuselage are computed for a quasi-static pull-up maneuver, a 
landing maneuver, and travel over runway bumps. These three load cases occur at user- 
specified fractions of gross take-off weight. For pitch control there is an option to use 
either elevators mounted on the horizontal tail (the conventional configuration) or elevons 
mounted on the trailing edges of the wing. The envelope of maximum bending moments is 
computed for all three load cases and is then used to determine the net stress resultants at 
each fuselage station. 

After the net stress resultants are determined at each fuselage station, a search is 
conducted at each station to determine the amount of structural material required to 
preclude failure in the most critical condition at the most critical point on the shell 
circumference. This critical point is assumed to be the outermost fiber at each station. 
Failure modes considered are tensile yield, compressive yield, local buckling, and gross 
buckling of the entire structure. A minimum gage restriction is also imposed as a final 
criteria. It is assumed that the material near the neutral fiber of the fuselage (with respect 
to longitudinal bending loads) is sufficient to resist the shear and torsion loads transmitted 
through the fuselage. For the shear loads this is a good approximation as the fibers 


farthest from the neutral axis will carry no shear, also, for beams with large fineness ratios 
(fuselage length/maximum diameter) bending becomes the predominant failure mode. 

The maximum stress failure theory is used for predicting yield failures. Buckling 
calculations assume stiffened shells behave as wide columns and sandwich shells behave as 
cylinders. The frames required for the stiffened shells are sized by the Shanley criterion. 
This criterion is based on the premise that, to a first-order approximation, the frames act 
as elastic supports for the wide column. 

There are a variety of structural geometries available for the fuselage. There is a 
simply-stiffened shell concept using longitudinal frames. There are three concepts with Z- 
stiffened shells and longitudinal frames; one with structural material proportioned to give 
minimum weight in buckling, one with buckling efficiency compromised to give lighter 
weight in minimum gage, and one a buckling-pressure compromise. Similarly, there are 
three truss-core sandwich designs, two for minima! weight in buckling with and without 
frames, and one a buckling-minimum gage compromise. 

It is assumed that the structural materials exhibit elasto-plastic behavior. Further, 
to account for the effects of creep, fatigue, stress-corrosion, thermal cycling and thermal 
stresses, options are available to scale the material properties of strength and Young’s 
modulus of elasticity. In the numerical results of this study, all materials were considered 
elastic and the full room-temperature material properties were used. 

Composite materials can be modeled with PDCYL by assuming them to consist of 
orthotropic lamina formed into quasi-isotropic (two-dimensionally, or planar, isotropic) 
laminates. Each of the lamina is assumed to be composed of filaments placed 
unidirectionally in a matrix material. Such a laminate has been found to give very nearly 
minimum weight for typical aircraft structures. 

Wing 

The wing structure is a multi-web box beam designed by spanwise bending and 
shear. The wing-fuselage carrythrough structure, defined by the wing-fuselage 
intersection, carries the spanwise bending, shear, and torsion loads introduced by the 
outboard portion of the wing 

The load case used for the wing weight analysis is the quasi-static pull-up 
maneuver. The applied loads to the wing include the distributed lift and inertia forces, and 
the point loads of landing gear and propulsion, if placed on the wing. Fuel may also be 
stored in the wing, which will relieve bending loads during the pull-up maneuver. 

The wing weight analysis proceeds in a similar fashion to that of the fuselage. The 
weight of the structural box is determined by calculating the minimum amount of material 
required to satisfy static buckling and strength requirements at a series of spanwise 
stations. The covers of the multi-web box are sized by buckling due to local instability 
and the webs by flexure-induced crushing. Required shear material is computed 
independently of buckling material. Aeroelastic effects are not accounted for directly, 
although an approximation of the magnitude of the tip deflection during the pull-up 


maneuver is made. For the carrythrough structure, buckling, shear, and torsion material 
are computed independently and summed. 

As for the fuselage, there are a variety of structural geometries available. There 
are a total of six structural concepts, three with unstiffened covers and three with truss- 
stiffened covers. Both cover configurations use webs that are either Z-stiffened, 
unflanged, or trusses. 


Geometry 

Fuselage 

The fuselage is assumed to be composed of a nose section, an optional cylindrical 
midsection, and a tail section. The gross density and fineness ratio are defined as 
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where W T o is the gross take-off weight, V B is the total fuselage volume, l B is the fuselage 
length, and D is the maximum fuselage diameter. The fuselage outline is defined by two 
power-law bodies of revolution placed back-to-back, with an optional cylindrical 
midsection between them (Figure 1 .) (For the present study, all nine transports used for 
validation of the analysis used the optional cylindrical midsection.) 

With the cylindrical midsection, integration gives the fuselage volume, fuselage 
planform area, and fuselage surface area as 
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respectively, where l t and h are the respective lengths to the start and end of the 
cylindrical midsection, and Pi and P2 are the respective powers that describe the nose and 
tail sections. Pi and P 2, again for the case of the cylindrical midsection, are found by 
solving the power-law equations for the volumes of the nose and tail sections, which are 
input from ACSYNT. The solution of these equations gives the respective nose and tail 
powers as 
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where V/ and V 2 are the corresponding nose and tail volumes. 

The horizontal tail is placed according to its quarter chord location as a fraction of 
the fuselage length. The distance from the nose to the tail is 


X HT ~ ^B^HT 


( 8 ) 


where Rht is the ratio of horizontal tail station to fuselage length. 

Propulsion may be either mounted on the fuselage or placed on the wing. In the 
case of fuselage mounted propulsion, the starting and ending positions of the propulsion 
unit are again calculated from their respective fractions of fuselage length as 
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where R P , and R P2 are the corresponding ratios of lengths to the leading and trailing edges 
of the fuselage engine pod to fuselage length. 

Similarly, the nose landing gear is placed on the fuselage as a fraction of vehicle 
length; the main gear, on the other hand, may be placed either on the fuselage as a single 
unit, also as a fraction of fuselage length, or on the wing in multiple units as will be 
described below. The positions of the respective nose and optional fuselage-mounted 
main gear are 
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where Rkg and Rkig are the corresponding length ratios for the nose gear and main gear 
stations to vehicle length. 

Wing 

The lifting planforms are assumed to be tapered, swept wings with straight leading 
and trailing edges. 

The wing loading is defined as 
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where Sp is the wing planform area. 

The wing is placed on the fuselage according to the location of the leading edge of 
its root chord, determined as a fraction of the fuselage length. The distance from the nose 
to the leading edge of the wing is 
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where Rle is the ratio of leading edge station to fuselage length. 

The first step in computing the wing weight is the determination of the geometry 
of the structural wing box. In terms of the input parameters Wpo, (W/S p ), aspect ratio 
( AR ), taper ratio ( Rtap ), and leading edge sweep (Ale), the dependent parameters wing 
area, span, root chord, tip chord, and trailing edge wing sweep are computed from 
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(Figure 2). It is assumed that specified portions of the streamwise (aerodynamic) chord 
are required for controls and high lift devices, leaving the remainder for the structural 
wing box. The portions of the leading and trailing edges that are left for non-structural 
use are specified as respective fractions C s> and C 5j of the streamwise chord. 

Determination of these chord fractions is accomplished through visual inspection of the 
wing planform. Measured at the theoretical root chord, the dimensions for the leading and 
trailing edges are 
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respectively. The intersection of this structural box with the fuselage contours determines 
the location of the rectangular carrythrough structure. The width of the carrythrough 
structure, w c , is defined by the corresponding fuselage diameter. 

The dimensions of the structural box and of the carrythrough structure are now 
determined (Figure 3). The structural semispan, b s , is assumed to lie on the quarter-chord 
line,_y, whose sweep is given by 
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Thus, 
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The streamwise chord at any point on the wing is given by 
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where C, is measured perpendicular to the vehicle longitudinal axis from the vehicle 
centerline toward the wingtip. Thus, the streamwise chord is the dimension of the wing 
parallel to the vehicle longitudinal axis. In particular, at the wing-fuselage intersection. 
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The structural root and tip chords are 
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respectively. In terms ofy, measured along the quarter chord from the wing-fuselage 
intersection toward the wingtip, the structural and total chords are given by 
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where the structural chord is defined as the dimension of the rectangular wing box 
measured parallel to the vehicle longitudinal axis. Computation of the widths of the wing 
box and total wing structure, as shown in Figure 3, is relatively complicated due to the 
geometry at the wingtip and the wing-fuselage intersection. For the portion of the wing 
between the wingtip and the wing-fuselage intersection, the respective widths of the wing 
box and total wing structure at any spanwise station^ are 


Z s (y) = r s cos(A s ) (30) 

Z(y) = r cos(A 5 ) (31) 


where Zs(y) and Z{y) are dimensions perpendicular to the structural semispan. 

The thickness of the wing box at any spanwise station^ is 

frR/y), 0 < y < b s (box structure) 

|r/^(0), y <0 (carrythrough structure) 

where R,(y) is determined as a linear interpolation between the root and tip thickness 
ratios. 

For the transports in the present study, all the fuel is carried within the wing 
structure. An option is also available to carry the fuel entirely within the fuselage, 
negating any bending relief in the wing (The high altitude drone, described in Appendix 
C, was modeled with a fuselage fuel tank.) The volume of the rectangular wing box 
structure is found as follows: 
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This equation is based on flat upper and lower surfaces and neglects the volume taken up 
by the structure. 


Loads 


Fuselage 

Fuselage loading is determined on a station-by-station basis along the length of the 
vehicle. Three types of fuselage loads are considered - longitudinal acceleration, tank 
pressure, and bending moment. In the present study, all three load types are assumed to 
occur simultaneously to determine maximum compressive and tensile loads at the outer 
shell fibers at each station. 

Bending loads applied to the vehicle fuselage are obtained by simulating vehicle 
pitch-plane motion during a quasi-static pull-up maneuver; a landing, and movement over 
a runway bump. Simplified vehicle loading models are used where it is assumed that: (1) 
fuselage lift forces (nominally zero for subsonic transports) are distributed uniformly over 
the fuselage plan area, (2) wing loading, determined independently, is transferred by a 
couple of vertical force and torque through the wing carrythrough structure, (3) fuselage 
weight is distributed uniformly over fuselage volume, (4) control surface forces and 
landing gear reactions are point loads, and (5) the propulsion system weight, if mounted 
on the fuselage, is uniformly distributed. A factor of safety (nominally 1 .5) is applied to 
each load case. The aircraft weight for each case is selected as a fraction of gross take-off 
weight. The resulting one-dimensional loading model is shown in Figure 4. All fuselage 
lift forces are assumed to be linear functions of angle of attack. Longitudinal bending 
moments are computed for each of the three loading cases and the envelope of the 
maximum values taken as the design loading condition. The bending moment computation 
is given in detail in Reference 4 and will only be summarized here. 

Considering first the pull-up maneuver loading, the motion is assumed to be a 
quasi-static pitch-plane pull-up of given normal load factor n (nominally 2.5 for transport 
aircraft). The vehicle is trimmed with the appropriate control surface (a horizontal tail for 
all nine transport used for validation in the present study), after which the angle of attack 
is calculated. 

Landing loads are developed as the aircraft descends at a given vertical speed, Vs, 
after which it impacts the ground; thereafter the main and nose landing gear is assumed to 
exert a constant, or optionally a ( 1 -cos(cot)). force during its stroke, Slc, until the aircraft 
comes to rest. The vehicle weight is set equal to the nominal landing weight. Wing lift as 
a fraction of landing weight is specified, which reduces the effective load the landing gear, 
carries. Likewise, the portion of total vehicle load the main gear carries is specified. No 
pitch-plane motion is considered during the landing. 

Runway bump loads are handled by inputting the bump load factor into the landing 
gear. This simulates the vehicle running over a bump during taxi. In a similar fashion to 
the landing, the wing lift as a fraction of gross take-off weight is specified, as is the 
portion of effective load input through the main gear. No pitch-plane motion is 
considered during the bump. 


Wing 

For the wing, only a quasi-static pull-up maneuver condition at load factor n is 
considered for determining loads. At each spanwise station along the quarter chord, from 
the wingtip to the wing-fuselage intersection , the lift load, center of pressure, inertia load, 
center of gravity, shear force and bending moment are computed. For the inertia load, it is 
assumed that the fuel weight Wft is distributed uniformly with respect to the wing volume 
so that the inertial load at y is (WFr/Vw)*V(y), where V(y) is the volume outboard ofy; this 
volume has centroid C g (y) with respect to station y. An estimate of the wing structural 
weight is included in Wrr for this calculation but the calculation is not redone when the 
actual structural weight has been computed. 

There is an option for either a trapezoidal or a Schrenk (Reference 16) lift load 
distribution along the wingspan; the trapezoidal distribution represents a uniform lift over 
the wing area (which has a trapezoidal planform) while the Schrenk distribution is an 
average of the trapezoidal distribution with an elliptical distribution, where the lift is zero 
at the wingtip and maximum at the wing-fuselage intersection. Prandtl has shown that a 
true elliptical lift load distribution will have a minimum induced drag, but a combination of 
the elliptical and trapezoidal distributions will give a better representation of actual aircraft 
loading (Reference 16) 

Plots of trapezoidal and Schrenk lift load distributions are shown in Figure 5. For 
the trapezoidal lift load distribution the lift load at y is (W/S)A T rap (y), where A T Mp(y) is 
the area outboard ofy; the centroid of this area is denoted C P {y) , where y is measured 

along the quarter chord. For the elliptical lift load distribution, the lift load matches the 
contour of an ellipse with the end of its major axis on the tip and the end of its minor axis 
directly above the wing-fuselage intersection. The area enclosed by the quadrant of the 
ellipse is set equal to the exposed area of the trapezoidal wing panel: 
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Thus the value of lift at^, Lell , the area of ellipse outboard ofy, Aell , and the center of 
pressure of lift outboard o fy, C Piu , for^ measured along the structural box may be 

determined as 
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respectively. 

For the Schrenk lift load distribution, the average of ArRAp(y) and A rn (v) is used to 
represent the composite area, while the average of C Pm Jy) and C P (y) is used to 

represent the composite center of pressure. 

Using the appropriate outboard area A(y) and center of pressure Cp(y), the shear 

force is 
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where ii e and n lg are the number of engines and landing gear mounted on the semispan, 
respectively; W t and are the weights of the i' h engine an i ,h landing gear, respectively, 

y t and _y lg are the locations of the i' h engine and i ,h landing gear, respectively; and 
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The bending moment is 
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Structural Analysis 
Fuselage 

Weight estimating relationships are now developed for the load-carrying fuselage 
structure. In addition, the volume taken up by the fuselage structure is also determined. 

Considering first the circular shell, the stress resultants in the axial direction caused 
by longitudinal bending, axial acceleration, and pressure at a fuselage station x are 
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respectively, where r=D 2 is the fuselage radius, A=nr* is the fuselage cross-sectional area, 
and P-2nr is the fuselage perimeter. In Equation 42, ly'-nr 3 is the moment of inertia of 
the shell divided by the shell thickness. In equation 43, for the case of fuselage-mounted 
propulsion, W s is the portion of vehicle weight ahead of station x if x is ahead of the inlet 
entrance, or the portion of vehicle weight behind x if x is behind the nozzle exit. In 
equation 44, P g is the limit gage pressure differential for the passenger compartment 
during cruise. The total tension stress resultant is then 
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if x is ahead of the nozzle exit, and 
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if ac is behind it. Similarly, the total compressive stress resultant is 
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if x is ahead of the inlet entrance, and 
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N x - = N Xt - 


if x is behind it. These relations are based on the premise that acceleration loads never 
decrease stress resultants, but pressure loads may relieve stress, if pressure stabilization is 
chosen as an option. The stress resultant in the hoop direction is 

N y =r P g K p (49) 


where K p accounts for the fact that not all of the shell material (for example, the core 
material in sandwich designs) is available for resisting hoop stress. 

The equivalent isotropic thicknesses of the shell are given by 
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for designs limited by compressive yield strength {Fey), ultimate tensile strength ( F „,), and 
minimum gage, respectively. In Equation 52, t mg is a specified minimum material thickness 
and K mg is a parameter relating t Sa to t mg which depends on the shell geometry. 

A fourth thickness that must be considered is that for buckling critical designs, t St , 

which will now be developed. The nominal vehicles of this study have integrally stiffened 
shells stabilized by ring frames. In the buckling analysis of these structures, the shell is 
analyzed as a wide column and the frames are sized by the Shanley criteria (Reference 9). 
Expressions are derived for the equivalent isotropic thickness of the shell required to 
preclude buckling, t Si , and for the smeared equivalent isotropic thickness of the ring 

frames required to preclude general instability, t F . The analysis will be restricted to the 
case of cylindrical shells. The major assumptions are that the structural shell behaves as an 
Euler beam and that all structural materials behave elastically. 

For the stiffened shell with frames concept, the common procedure of assuming 
the shell to be a wide column is adopted. If the frame spacing is defined as d and Young’s 
modulus of the shell material is defined as E, the buckling equation is then 
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or, solving for i St 
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Values of the shell efficiency e for the various structural concepts are given in 
Table 2; the structural shell geometries available are simply stiffened, Z-stiffened, and 
truss-core sandwich. We next size the frames to prevent general instability failure. The 
Shanley criterion is based on the premise that the frames act as elastic supports for the 
wide column, this criterion gives the smeared equivalent thickness of the frames as 
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where C/ris Shanley’s constant, K Ft is a frame geometry parameter, and Ef is Young’s 
modulus for the frame material (See Reference 3 for a discussion of the applicability of 
this criterion and for a detailed derivation of the equations presented here.) If the 
structure is buckling critical, the total thickness is 
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Minimizing t with respect to d results in 
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where Pf is the density of the frame material and p is the density of the shell material, so 
that the shell is three times as heavy as the frames. 

Frame-less sandwich shell concepts may also be used. For these concepts, it is 
assumed that the elliptical shell buckles at the load determined by the maximum 
compressive stress resultant N x on the cylinder. The buckling equation for this frame-less 
sandwich shell concepts is 
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where m is the buckling equation exponent. Or, solving for t Sf 
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This equation is based on small deflection theory, which seems reasonable for sandwich 
cylindrical shells, although it is known to be inaccurate for monocoque cylinders. Values 
of m and e may be found, for example in References 5 and 6 for many shell geometries. 
Table 2 gives values for sandwich structural concepts available in PDCYL, numbers 8 and 
9, both of which are truss-core sandwich. The quantities N x \ r, and consequently , will 
vary with fuselage station dimension x. 

At each fuselage station x, the shell must satisfy all failure criteria and meet all 
geometric constraints. Thus, the shell thickness is selected according to compression, 
tension, minimum gage and buckling criteria, or 
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If t s = t s , the structure is buckling critical and the equivalent isotropic thickness of the 
frames, t F , is computed from Equation 59. If t s > , the structure is not buckling critical 

at the optimum frame sizing and the frames are re-sized to make t s = i s> . Specifically, a 
new frame spacing is computed from Equation 54 as 
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and this value is used in Equation 55 to determine t F . 

The total thickness of the fuselage structure is then given by the summation of the 
smeared weights of the shell and the frames 
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The shell gage thickness may be computed from t g =t s IK mg . The ideal fuselage structural 
weight is obtained by summation over the vehicle length 

W, = 27tl(p t Sf + p F i F )r l Ax l (66) 


where the quantities subscripted i depend on x 

We next discuss the derivation of the structural geometry parameters shown in 
Table 2. The Z-stiffened shell, typical of modem transport aircraft, will be used as an 
example of skin-stringer-frame construction. Using Reference 5 and Figure 6, the 
equivalent isotropic thickness of the smeared skin and stringers is 
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Since only the skin is available for resisting pressure loads, 
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For minimum gage designs, if t s >t w then and 
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On the other hand, if l s <t w then t s =t mg and 
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Equations 68, 70, and 72 show that for both pressure loading critical and minimum gage 
limited structure, ( bjb s ) and (ljt s ) should be as small as possible (i.e. no stringers). As an 
option in PDCYL, all of the detailed shell dimensions shown in Figure 6 are computed and 
output at each fuselage station. 

In practice, a typical design will be influenced by bending and pressure loads and 
by the minimum gage constraint, and thus a compromise is necessary. If buckling is of 
paramount importance, then a good choice is ( bjb s )=0.$l and (tjt s )- 1.06 because this 
gives the maximum buckling efficiency for this concept, namely e=0.91 1 (Reference 5). 
From Equations 68 and 72, 


^ = ^ = 1 + (1.6X0.8 7 )( 1 .06) = 2.475 


(73) 


This is concept 3 in Tables 1 and 2. If pressure dominates the loading condition, then 
(bjb t )= 0.6 and (tjt s )= 0.6 is a reasonable choice, giving e=0.76, AT p = 1 .576, and 
AT m g=2.628; this is concept 5. For minimum gage dominated structure, the geometry 
{bjb t y*0.5% and (tjt s )=0.90 gives concept 6. 


The geometry of the truss-core sandwich shell concept is shown in Figure 7. The 
equivalent isotropic shell thickness of this concept is 


h ~ 


2 +- 


1 


t, cos(vy)J 


(74) 


Reference 5 shows that the optimum buckling efficiency is obtained for {tjtjjr 0.65 and 
V=55°. This gives £=0.4423, K mg =4.%20, and K p - 3. 132, concept 8 in Tables 1 and 2. To 
get a design that is lighter for minimum gage dominant structure, a geometry is chosen 
that places equal thickness material in the face sheets and the core; the choice of 
(JtJtfjr 10 and vj/=45° gives structural concept 9. These calculations assume that the face 
sheets and core are composed of the same material and are subject to the same minimum 
gage constraint. 

Since the preceding analysis gives only the ideal weight, Wj, the non-optimum 
weight, Wso, (including fasteners, cutouts, surface attachments, uniform gage penalties, 
manufacturing constraints, etc.) has yet to be determined. The method used will be 
explained in a later section. 

Wing 

Using the geometry and loads applied to the wing developed above, the structural 
dimensions and weight of the structural box may now be calculated. The wing structure is 
assumed to be a rectangular multi-web box beam with the webs running in the direction of 
the structural semispan Reference 5 indicates that the critical instability mode for 
miltiweb box beams is simultaneous buckling of the covers due to local instability and of 
the webs due to flexure induced crushing. This reference gives the solidity (ratio of 
volume of structural material to total wing box volume) of the least weight multi-web box 
beams as 


I = e 


f M ' 
\Z s t'Ej 


(75) 


where e and e depend on the cover and web geometries (Table 3), Mis the applied 
moment, t is the thickness, E is the elastic modulus, and Zs is obtained from Reference 5. 
The solidity is therefore 


WeudO') 

pZ s t 


( 76 ) 


where Wbend is the weight of bending material per unit span and p is the material density. 
W bend is computed from Equations 75 and 76. The weight per unit span of the shear 
material is 




(77) 


where F s is the applied shear load and Os is the allowable shear stress. The optimum web 
spacing (Figure 8) is computed from (Reference 2) 

J*c_ 

4'c-3 

(78) 


d w =t 


(l-2e c ) f M 
(\-e c )j^X Z s rE 


where subscripts W and C refer to webs and covers, respectively. The equivalent isotropic 
thickness of the covers and webs are 



respectively, and the gage thicknesses are 

(81) 
(82) 



Values of e, e, Ec, E c , £w, K gr , and K gc , are found in Table 3 for various structural 

concepts (Reference 5). If the wing structural semispan is divided into N equal length 
segments, the total ideal weight is the wing box structure is 


W, 


BOX 


_ 2*s y ( 

N 


W' 

"BEH'D, 


+ W' ^ 

T "shear, f 


(83) 


The wing carrythrough structure consists of torsion material in addition to bending 
and shear material The torsion material is required to resist the twist induced due to the 
sweep of the wing. The bending material is computed in a similar manner as that of the 
box except that only the longitudinal component of bending moment contributes. Letting 
t(ni(y=-0) and 0 ) , 


2c 


r cos(A 5 ) Y 

l J 


(84) 


The weight of the bending material is then 


^S£N'D r ~ P^CQ«'0»C 


(85) 


where vt > c is the width of the carrythrough structure. (When the wing-fijselage intersection 
occurs entirely within the cylindrical midsection, as is the case with all nine transport used 
for validation in the present study, wc=D.) The quantities dw, tw, and t c , are computed in 
the same manner as for the box. The weight of the shear material is 


W Q 


SHEARr 



( 86 ) 


where F s =F s {0). 

The torque on the carrythrough structure is 


T = M 0 sin(A 5 ) 


( 87 ) 


and the weight of the torsion material is then 


W - 

TORSION c 


P7X'q Qg? ) W C 


( 88 ) 




S 


Finally, the ideal weight of the carrythrough structure is computed from a summation of 
the bending shear and torsion material, or 

W c = W BEND c + ^ SHEAR c + ^TORSION c (89) 


As in the case of the fuselage structural weight, non-optimum weight must be 
added to the ideal weight to obtain the true wing structural weight. The method used will 
be discussed below. 

The static deflection of the wingtip under the pull-up maneuver is also determined. 
Using the moment-area method applied to an Euler beam (Reference 14), the deviation of 
point B on the deflected surface from the tangent drawn from another point A on the 
surface is equal to the area under the M/(E1) diagram between A and B multiplied by the 
distance to the centroid of this area from B, 

(90) 


where 0 is the angular displacement of the beam and y is the longitudinal axis of the beam. 
For the case of a wing with trapezoidal planform, the longitudinal axis,^, will lie along the 
quarter-chord line (Figure 3). For a wing with a horizontal unloaded configuration, the 
tangential deviation, 1 B a, will equal the true vertical tip displacement (assumed to be the 
case). Only the wing cover contributes to the bending resistance, while the webs offer 
similar shear stiffness. The wing area moment of inertia, /, at any structural semispan 
station^ is determined with the Parallel Axis theorem, as cover thickness is small when 
compared with total wing thickness. 


Regression Analysis 
Overview 

Using fuselage and wing weight statements of nine subsonic transports, a relation 
between the calculated load-bearing structure weights obtained through PDCYL and the 
actual load-bearing structure weights, primary structure weights, and total weights is 
determined using statistical analysis techniques. A basic application which is first 
described is linear regression, wherein the estimated weights of the aircraft are related to 
the weights calculated by PDCYL with a straight line, y=mx+b, where y is the value of the 
estimated weight, m is the slope of the line, x is the value obtained through PDCYL, and b 
is the ^-intercept. This line is termed a regression line, and is found by using the method 
of least squares, in which the sum of the squares of the residual errors between actual data 
points and the corresponding points on the regression line is minimized. Effectively, a 
straight line is drawn through a set of ordered pairs of data (in this case nine weights 
obtained through PDCYL and the corresponding actual weights) so that the aggregate 
deviation of the actual weights above or below this line is minimized. The estimated 
weight is therefore dependent upon the independent PDCYL weight. 

As an example, if the form of the regression equation is linear, the estimated 
weight is 


y t s<= mx c*,c +b 


( 91 ) 


where m is the slope, b is the intercept, and x ca ic is the weight PDCYL calculates. The 
resulting residual to be minimized is 


£ = 2 > 

1 = 1 


actual 



( 92 ) 


or 


£ = IO« tual, ~ m X calc, - bf 
i = l 


( 93 ) 


where y actual is the actual component weight and n is the number of aircraft whose data is 
to be used in the fit. By taking partial derivatives of the residual error with respect to both 
m and b, equations for the values of these two unknown variables are found to be 


m = 


(94) 


n n n 

»Z - Z Z 

r = l i=J » = ) 
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»Z x Lc.--(Z *«*,)’ 

1 = 1 i=l 


6 = y^ - w x^ x,y> = mean values of x andy 


(95) 


Of key importance is the degree of accuracy to which the prediction techniques are 
able to estimate actual aircraft weight. A measure of this accuracy, the correlation 
coefficient, denoted R, represents the reduction in residual error due to the regression 
technique. R is defined as 


R = 



( 96 ) 


where £, and E r refer to the residual errors associated with the regression before and 
after analysis is performed, respectively. A value of R= 1 denotes a perfect fit of the data 
with the regression line Conversely, a value of R = 0 denotes no improvement in the data 
fit due to regression analysis. 

There are two basic forms of equations which are implemented in this study. The 
first is of the form 


y^= mx caic 


( 97 ) 


The second general form is 


y, s , = m x °aic + b W 

The first form is a simplified version of the linear example as discussed above, with the y- 
intercept term set to zero. However, because the second general equation is not linear, 
nor can it be transformed to a linear equation, an alternative method must be employed. 

In order to formulate the resulting power-intercept regression equation, an iterative 
approach developed by D. W. Marquardt is utilized (Reference 15). This algorithm starts 
at a certain point in space, and, by applying the method of steepest descent, a gradient is 
obtained which indicates the direction in which the most rapid decrease in the residual 
errors will occur. In addition, the Taylor Series method produces a second similar vector. 
Interpolation between these two vectors yields a direction in which to move the point in 
order to minimize the associated error. After several iterations, the process converges to a 
minimum value. It should be noted that there may be several local minimums and there is 
no guarantee that the method converges to the global one. 


Fuselage 

The analysis above is used to develop a relationship between weight calculated by 
PDCYL and actual wing and fuselage weights. The data was obtained from detailed 
weight break-downs of nine transport aircraft (References 12, 13, 17-19) and is shown in 
Table 4 for the fuselage. Because the theory used in the PDCYL analysis only predicts the 
load-carrying structure of the aircraft components, a correlation between the predicted 
weight and the actual load-carrying structural weight, primary weight, as well as the total 
weight of the fuselage was made. 

Structural weight consists of all load-carrying members including bulkheads and 
frames, minor frames, covering, covering stiffeners, and longerons. For the linear curve- 
fit, the resulting regression equation is 


^actual = 1-3313 *=0.99624 


(99) 


This shows that the non optimum factor for fuselage structure is 1.3313; in other words, 
the calculated weight must be increased by about 33% to get the actual structural weight. 
For the alternative power-intercept curve fitting analysis, the resulting load-carrying 
regression equation is 


= 4019.9 + 0.0 1 952 W^ 073 *=0.99624 


( 100 ) 


To use either of these equations to estimate total fuselage weight, non-structural weight 
items must be estimated independently and added to the structural weight. 

Primary weight consists of all load-carrying members as well as any secondary 
structural items such as joints fasteners, keel beam, fail-safe straps, flooring, flooring 
structural supplies, and pressure web. It also includes the lavatory structure, galley 
support, partitions, shear ties, tie rods, structural firewall, torque boxes, and attachment 
fittings. The linear curve fit for this weight yields the following primary regression 
equation 


^= 18493^ *=0.98471 


( 101 ) 


The primary power-intercept regression equation is 


=553 1.9 + 0.01 9094 *= 0.98934 


( 102 ) 


The total fuselage weight accounts for all members of the body, including the 
structural weight and primary weight. It does not include passenger accommodations, 
such as seats, lavatories, kitchens, stowage, and lighting; the electrical system; flight and 
navigation systems; alighting gear; fuel and propulsion systems; hydraulic and pneumatic 



systems; the communication system; cargo accommodations, flight deck 
accommodations; air conditioning equipment; the auxiliary power system; and 
emergency systems. Linear regression results in the following total fuselage weight 
equation 


*2-5178^ £=098549 


( 103 ) 


The total fuselage weight power-intercept regression equation is 


W xntal = 6877.8 + 0.099541 £=0.9878 


(104) 


Plots of actual fuselage component weight verses PDCYL-calculated weight, as 
well as the corresponding linear regressions, is shown in Figures 9-11. 

Wing 

The same analysis was performed on the wing weight for the sample aircraft and is 
shown in Table 5. The wing box, or load-carrying structure, consists of spar caps, 
interspar coverings, spanwise stiffeners, spar webs, spar stiffeners, and interspar ribs. The 
wing box linear regression equation is 


^ = 1-0379 W calc R= 0.98482 


(105) 


so that the nonoptimum factor, is 1 .0379. Power-intercept regression results in 


= -930.36 + 1.13 1 W*™ 12 £= 0.98537 


(106) 


Wing primary structural weight includes all wing box items in addition to auxiliary 
spar caps and spar webs, joints and fasteners, landing gear support beam, leading and 
trailing edges, tips, structural firewall, bulkheads, jacket fittings, terminal fittings, and 
attachments. Linear regression results in 


^ = 142171^ £=0.99285 


(107) 


Power-intercept regression yields 


=-908.14 + 1.8928 W?™ 


£= 0.9929 


(108) 


The total wing weight includes wing box and primary weight items in addition to 
high-lift devices, control surfaces and access items. It does not include the propulsion 
system, fuel system, and thrust reversers; the electrical system; alighting gear; hydraulic 
and pneumatic systems; anti-icing devices; and emergency systems. The resulting total 
weight linear regression equation is 


= 1.83530^ *=0.98857 (109) 

The power-intercept equation for total wing weight is 

= 657.33 + 1.662 4 *=0.98862 (110) 


Plots of actual wing component weight verses PDCYL-calculated weight, as well 
as the corresponding linear regressions, is shown in Figures 12-14. 

Discussion 

Both fuselage and wing weight linear regressions give excellent correlation with 
the respective weights of existing aircraft, as evidenced by the high values of the 
correlation coefficient, *. 

There are, however, two results of the correlations which are somewhat puzzling. 
First, the constant term in two of the wing power-intercept equations is negative; this will 
usually numerically insignificant Second, the exponents in the fuselage power-intercept 
equations are significantly different than one; this may be due to convergence to a local, 
rather than the global, minimum. 

Because estimates of non-load-bearing primary structure are generally not 
available at the conceptual design stage, and because non-primary structure is probably 
not well-estimated by a non-optimum factor, Equation 101 and Equation 107 are 
recommended for estimating the primary structural weights of the respective transport 
fuselage and wing structures (Figure 10, Figure 13). 


REFERENCES 


1. Ardema, MD., Body Weight of Hypersonic Aircraft: Parti. NASA TM 101028, 
Oct. 1988. 


2. Ardema, M. D. and Williams, L. J., Transonic Transport Study - Structures and 
Aerodynamics, NASA TM X-62,157, May 1972. 


3. Ardema, M. D., Structural Weight Analysis of Hypersonic Aircraft, NASA 
TND-6692, Mar. 1972. 


4. Ardema, M. D., Analysis of Bending Loads of Hypersonic Aircraft, NASA 
TMX-2092, 1970. 


5. Crawford, R. F. and Bums, A B , Strength, Efficiency, and Design Data for 
Beryllium Structures, ASD-TR-6I-692, Feb. 1962. 


6. Crawford, R. F. and Bums, A. B., Minimum Weight Potentials for Stiffened Plates 
and Shells, AIAA Journal, vol. 1, no. 4, Apr. 1963, pp. 879-886. 

7. Vinson, J. R. and Shore, S., Minimum Weight Web Core Sandwich Panels Subjected 
to Uniaxial Compression, Journal of Aircraft, vol. 8, no. 11, Nov. 1971, pp. 843-847. 

8. Vinson, J. R., Optimum Design of Composite Hex-Cell and Square Cell Honeycomb 
Sandwich Panels Subjected to Uniaxial Compression, AIAA Journal, vol. 24, no. 10, 
Oct. 1987, pp. 1690-1696. 

9. Shanley, F. R., Weight-Strength Analysis of Aircraft Structures, Second Edition, 
Dover, New York, 1960. 


10. Timoshenko, S. P. and Young, D. H., Theory of Structures, McGraw-Hill, 1945. 

11. Moore, M. and Samuels, J., ACSYNT Aircraft Synthesis Program - User’s Manual, 
Systems Analysis Branch, NASA- Ames Research Center, Sep. 1990. 


12. Niu, M. C.-Y., Airframe Structural Design: Practical Design Information and Data 
on Aircraft Structures , Conmilit Press, 1991. 

13. York, P. and Labell, R. W., Aircraft Wing Weight Build-Up Methodology with 
Modification for Materials and Construction Techniques, NASA CR 166173, Sep. 
1980. 


14. Popov, E P., Mechanics of Materials, Second Edition, Prentice-Hall, New Jersey, 
1976. 


15. Marquardt, D W., Least Squares Estimation of Nonlinear Parameters, a Computer 
Program in FORTRAN IV Language, IBM SHARE Library, Distribution Number 
309401, August 1966. 


16 McCormick, B W., Aerodynamics, Aeronautics, and Flight Mechanics, John Wiley 
& Sons, 1979. 


17. The Boeing Company, Commercial Airplane Division, Weight Research Group, 
Thomas, R. B. and Parsons, S. P., Weight Data Base, Doc. No. D6-23204 TN, 1968. 

18. McDonnell Douglas Aircraft Company, Detailed weight statement for MD-1 1 
transport aircraft, June 1987. 

19. McDonnell Douglas Aircraft Company, Detailed weight statement for MD-80 
transport aircraft, July 1991. 

20. Megson, T. H. G., Aircraft Structures for Engineering Students, Second Edition, 
Halsted Press, 1990. 


TABLE 1. FUSELAGE STRUCTURAL GEOMETRY CONCEPTS. 


(KCON Sets Concept Number) 

2. Simply Stiffened Shell, Frames, Sized for Minimum Weight in Buckling 

3. Z-Stiffened Shell, Frames, Best Buckling 

4. Z-Stiffened Shell, Frames, Buckling-Minimum Gage Compromise 

5. Z-Stiffened Shell, Frames, Buckling-Pressure Compromise 

6. Truss-Core Sandwich, Frames, Best Buckling 

8. Truss-core Sandwich, No Frames, Best Buckling 

9. Truss-Core Sandwich, No Frames, Buckling-Minimum Gage -Pressure Compromise 


TABLE 2. FUSELAGE STRUCTURAL GEOMETRY PARAMETERS. 


BSBaaBSSBil 

m 

8 


K, 

K* 

2 

2 

0.656 

2.463 

2.463 

0.0 

3 

2 

0.911 

2.475 

2.475 

0.0 

4 

2 

0.760 

2.039 

1.835 

0.0 

5 

2 

0.760 

2.628 

1.576 

0.0 

6 

2 

0.605 

4.310 

3.965 

0.459 

8 

1.667 

0.4423 

4.820 

3.132 

0.405 

9 

1.667 

0.3615 

3.413 

3.413 

0.320 















TABLE 3. WING STRUCTURAL COEFFICIENTS AND EXPONENTS. 


COVERS 

WEBS 

e 

e 

e 

e* 

Sw 



UNSTIFF 

TRUSS 

2.25 

0.556 

3.62 

3 

0.605 

1.000 

0.407 

UNSTIFF 

UNFLANGED 

2.21 

0.556 

3.62 

3 

0.656 

1.000 

0.505 

UNSTIFF 

Z-STIFF. 

2.05 

0.556 

3.62 

3 

0.911 

1.000 

0.405 

TRUSS 

TRUSS 

2.44 

0.600 

1.108 

2 

0.605 

0.546 

0.407 

TRUSS 

UNFLANGED 

2.40 

0.600 

1.108 

2 

0.656 

0.546 

0.505 

TRUSS 

Z-STIFF. 

2.25 

0.600 

1.108 

2 

0.911 

0.546 

0.405 

























TABLE 4. FUSELAGE WEIGHT BREAKDOWNS FOR NINE TRANSPORT AIRCRAFT. 



Weight, lb. 

Aircraft 

PDCYL 

Load-Carrying 

Structure 

Primary Structure 

Total Structure 

720 

6622 

9013 

13336 

19383 

727 

5848 

8790 

12424 

17586 

737 

3425 

5089 

7435 

11831 

747 

28507 

39936 

55207 

72659 

880 

8612 

8705 

9452 

13507 

DC-8 

9517 

13312 

18584 

24886 

MD-11 

20608 

25970 

34999 

54936 

MD-83 

7488 

9410 

11880 

16432 

L-1011 

21762 

28352 

41804 

52329 































TABLE 5. WING WEIGHT BREAKDOWNS FOR NINE TRANSPORT AIRCRAFT. 



Weight, lb. 

Aircraft 

PDCYL 

Load-Carrying 

Structure 

Primary Structure 

Total Structure 

720 

13588 

11747 

18914 

23528 

727 

8336 

8791 

12388 

17860 

737 

5706 

5414 

7671 

10687 

747 

49285 

50395 

68761 

88202 

880 

12319 

9241 

14144 

17585 

DC-8 

21892 

19130 

27924 

35330 

MD-11 

30561 

35157 

47614 

62985 

MD-83 

6888 

8720 

11553 

15839 

L-1011 

24234 

28355 

36101 

46233 
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Figure 5. Spanwise Wing Lift-load Distributions. 



Figure 6. Typical Z-Stiffened Shell Geometry. 



Figure 7. Truss-Core Sandwich Geometry 



Figure 8 Wing Structural Concept. 
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Figure 10. Fuselage Primary Structure and Linear Regression. 
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Figure 1 1 . Fuselage Total Structure and Linear Regression. 
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Figure 13. Wing Primary Structure and Linear Regression. 
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Figure 14. Wing Total Structure and Linear Regression. 


Page 1 



DDBody 

■ N-0 Fuselage 
^ Frames 
m Shell 








PDCYL Wing Weight, lb. 


STRAT07.XLC 



Box length as a fraction of root chord 
Figure Variation in wing weight as a function of structural box length. 
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APPENDIX A: User’s Manual, Example 


Description 

The purpose of this appendix is to give a detailed example of the input procedure 
used to allow PDCYL to calculate fuselage and wing weights for a sample transport 
aircraft during an ACSYNT run. A sample output from PDCYL will also be given. The 
Boeing 747-2 IP will be used for the example. The layout of the 747-2 IP is shown in 
Figure 15. The weights of the load-carrying portions of the fuselage and wing box for the 
747-2 IP will be calculated by PDCYL and scaled by the respective nonoptimum factors 
developed earlier to give estimates for the weights of the fuselage and wing. A 
comparison between methods currently used by ACSYNT to estimate fuselage and wing 
weights and PDCYL output will be made with the corresponding actual weights of the 
747-21P. 


Input 


PDCYL requires input from both the existing ACSYNT data structure and an 
additional namelist containing data required by PDCYL which are not contained within the 
current ACSYNT format. There are three steps to run PDCYL within ACSYNT. First, 
the aircraft type is specified in the ACSYNT Control input. Currently the Transport 
Aircraft type is used. Second, data within ACSYNT module namelists is required. The 
ACSYNT Geometry, Trajectory, and Weights modules supply data for PDCYL execution. 
PDCYL uses the WING, HTAIL, VTAIL, FUS, WPOD, and FPOD namelists from the 
Geometry module. From the Trajectory module, the TRDATA namelist is used. From 
the Weights module the OPTS namelist is used. Third, data from the PDCYLIN namelist 
is used. 

Variables used from ACSYNT namelists and the PDCYLIN namelist are given in 
Table 6 and Table 7, respectively. Default values for all variables are also given. These 
default values match the Boeing 747-2 IP. Key configuration parameters are given for 
each of the nine aircraft used in the validation study in Table 8. An example of the 
PDCYLIN namelist input for the 747-2 IP is shown in Figure 16. 

Output 

PDCYL weights output begins with the wing box and carrythrough structure 
analysis. The wing is sized during a quasi-static pull-up maneuver where the load factor is 
set equal to the ultimate load factor (nominally 3.75). Wing output contains three parts. 
First is the overall geometrical configuration. Second is a detailed station-by-station 
bending, shear, and torsion analysis and corresponding geometrical sizing along the span. 
Third is the detailed geometrical layout, loading, and weight breakdown of the 
carrythrough structure, weight breakdown of the wing components and deflection of the 
wingtip. This wing weight is multiplied by the non-optimum factor and returned to 


ACSYNT. An example of the PDCYL wing weight output for the 747-2 IP is shown in 
Figure 17. 

Next, the fuselage is analyzed. Fuselage output contains four parts. First is the 
overall geometrical layout and weight breakdown. Second is a station-by-station bending, 
shear and axial stress analysis. Up to three load cases are investigated In order they are a 
quasi-static pull-up maneuver, a landing maneuver, and travel over runway bumps. Third, 
the envelope of worst case loading is shown for each station, from which the shell and 
frames are sized. Corresponding unit weight breakdowns are also given. As an option, 
the detailed geometric configuration at each station may be output. Fourth, weights 
summaries are given for the top and bottom sections of the fuselage (nominally the same). 
These summaries are then averaged to give the weight summary of the entire fuselage. 

The fuselage weight, including the corresponding non-optimum factor, is returned to 
ACSYNT. An example of the PDCYL fuselage weight output for the 747-21 P is shown 
in Figure 18 

Figure 19a shows a comparison between fuselage weight estimates from the 
Sanders equation, the Air Force equation, and PDCYL with the actual fuselage weight of 
the 747-2 IP. Figure 19b shows a similar comparison for the wing weight. SLOPE and 
TECH factors were set to one for the comparisons in Figures 19a and 19b. 


TABLE 6. ACSYNT VARIABLES. 


1. Geometry module 


Namelist WING 


Variable 

T>pe 

Dinicnsio 

Description 

Units/Comment 

Default (747) 

SWEEP 

float 

1 

Sweep of wing. 

Degrees 

37.17 

KSWEEP 

integer 

1 

l-> Referenced to the leading edge. 

2 -> Referenced to the quarter chord. 
3-> Referenced to the trailing edge. 

2 

AR 

float 

1 

Aspect Ratio of wing. 


6.96 

TAPER 

float 

1 



0.2646 

TCROOT 

float 

1 

Thickness-to-chord ratio at the root 


0.1794 

TCTIP 

float 

1 

Thickness-to-chord ratio at the tip 


0.078 

ZROOT 

float 

1 

Elevation of MAC above fuselage reference 
plane, measured as a fraction of 
the local fuselage radius. 

-0.1 

AREA 

float 

1 

Planform Area of wing. 

ft A 2 

5469 

DIHED 

float 

1 


Degrees 

7 

XWING 

float 

1 

Ratio of distance measured from nose to 
leading edge of wing to total fuselage length. 

0.249 


Namelist HTAIL (horizontal tail) 


Variable 

Type 

Dimensio 

Description 

Units/Comment 

Default (747) 

SWEEP 

float 

1 

Sweep of tail 

Degrees 

34.29 

KSWEEP 

integer 

1 

l-> Referenced to the leading edge. 

2 -> Referenced to the quarter chord. 
3-> Referenced to the trailing edge. 

2 

AR 

float 

1 

Aspect Ratio of the horizontal wing. 

(span) A 2/Area 

3.625 

TAPER 

float 

1 

Taper Ratio of the horizontal wing. 


0.25 

TCROOT 

float 

1 

Thickness-to-chord ratio at the root 


0.11 

TCTIP 

float 

1 

Thickness-to-chord ratio at the tip. 


0.08 

ZROOT 

float 

1 

Elevation of MAC above fuselage re: 
plane, measured as a fraction of 
the local fuselage radius. 

erence 

0.69 

AREA 

float 

1 

Planform Area of the horizontal win ft A 2 

1470 

XHTAIL 

float 

1 

Position for trailing edge of tail root 
chord. If ZROOT LE. I, then XHTAIL is 
given as a fraction of body length. Else, 
XHTAIL is given as a fraction of the local 
vertical tail chord. 

1 


















































































Namelist VTAIL (vertical tail) 


Variable 

Type 

Dimensio 

Description Units/Comment 

Default (747) 

SWEEP 

float 

1 

Sweep of vertical tail Degrees 

45.73 

KSWEEP 

integer 

1 

1 -> Referenced to the leading edge. 

2 ->Referenced to the quarter chord. 
3-> Referenced t the trailing edge 

2 

AR 

float 

1 

Aspect Ratio of vertical tail. 

(span) A 2/Area 

1.247 

TAPER 

float 

1 

Taper Ratio of vertical tail. 

fntys fjsTy, k R WR 1 

0.34 

TCROOT 

float 

1 

Thickness-to-Chord ratio at root. 

0.1298 

TCTIP 

float 

1 

Thickness-to-Chord ratio at tip. 

0.089 

ZROOT 

float 

1 

Elevation of MAC above fuselage reference 
plane, measured as a fraction of 
the local fuselage radius. 

0.6 

AREA 

float 

1 

Planform Area of vertical tail. ft. A 2 

830 


Namelist FUS (fuselage) 


Variable 

Type 

Dimensio 

Description 

Units/Comment 

Default (747) 

FRN 

float 

1 

Fineness Ratio of the nose section. 

(Length/Diameter) 

2.13 

FRAB 

float 

1 

Fineness Ratio of after-body section 

(Length/Diameter) 

3.29 

BODL 

float 

1 

Length of fuselage 

ft. 

225.167 

BDMAX 

float 

1 

Maximum diameter of fuselage 

ft. 

20.2 


Namelist WPOD (wing-mounted propulsion pod) 


Variable 


Type Dimensio 


Description 


Units/Comment Default (747) 


DIAM 


float 


1 


Engine Diameter. 


ft. 


6.2 


LENGTH 


float 


Length of engine pod. 


ft. 


15 


float 


X location of nose of pod relative to leading 
edge of wing, given as a fraction of local 
chord of wing. (>0 if face of pod is behind 
leading edge of wing.) 


-0.631 


float 


Y location of center of pod, given as a 
fraction of semispan, measured from 
body centerline. 


0.241 


float 


Z location of center of pod above wing local 
chord, given as fraction of maximum 
pod diameter. 


-0.83 


SWFACT 


float 


Wetted Area multiplier. 


















































































Namelist FPOD (fusela 

ge-mounted propulsion pod) 

Variable 

Type 

Dimensio 

Description Units/Comment 

Default (747) 

DIAM 

float 

1 

Engine Diameter. 

ft. 

N/A 

LENGTH 

float 

1 

Length of engine pod. 

ft. 

N/A 

SOD 

float 

1 

Stand-off-distance, the distance from the 
pod wall to the fuselage wall, given as a 
fraction of maximum pod radius. 

N/A 

THETA 

float 

1 

Angular orientation of pod, THETA 
measured positive up from 
the horizontal reference plane. 

Degrees 

N/A 

X 

float 

1 

X location of nose relative to nose of 
pod, given as a fraction of body 
length. 

N/A 


2. Trajectory module 


Namelist TRDATA (used for load factors) 


Variable 

Type 

Dimensio 

Description 

Units/Comment 

Default (747) 

DESLF 

float 

1 

Design load factor 

N/A 

2.5 

ULTLF 

float 

1 

Ultimate load factor, usually 
1.5*DESLF 

N/A 

3.75 


3. Weights module 
Namelist OPTS 


Variable 


Type Dimensio 


Description 


Units/Comment Default (747) 


WGTO 


float 


1 


Gross-take-off-weight. 


lb. 


713000 


WE 


float 


Total Weight of propulsion system, 
(includes both wing and fuselage 
mounted engines.) 


Ib 


44290 






















































TABLE 7. PDCYL VARIABLES. 


Namelist PDCYLIN 


Wing 


Material Properties 


Variable 

Type 


Description 

Units/Comment 

Default (747) 

PS 

float 

1 

Plasticity Factor. 


1 

TMGW 

float 

1 

Min. gage thickness for the wing 

inches 

0.2 

EFFW 

float 

1 

Buckling efficiency of the web. 



0.656 

EFFC 

float 

1 

Buckling efficiency of the covers 


1.03 

ESW 

float 

1 

Young's Modulus for wing material 

psi 

1.07E+07 

FCSW 

float 

1 

Ult. compressive strength of wing 

psi 

54000 

DSW 

float 

1 

Density of the wing material. 

lb./in A 3 

0.101 

KDEW 

float 

1 

Knock-down-factor for Young's Modulus. 

1 

KDFW 

float 

1 

Knock-down-factor for Ultimate strength. 

1 

Geometric Parameters 

Variable 

Type 

Dimensio 

Description 

Units/Comment 

Default (747) 

ISTAMA 

integer 

1 

1 -> the position of the wing is unknown. 

2 -> the position of the wing is known 

2 

CS1 

float 

1 

Position of structural wing box from 
leading edge as % of root chord. 

0.088 

CS2 

float 

1 

Position of structural wing box from 
trailing edge as % of root chord. 

0.277 

Structural Concept 

Variable 

Type 

Dimensio 

Description 

Units/Comment 

Default (747) 

CLAQR 

float 

1 

Ratio of body lift to wing lift 

For subsonic aircraf 
CLAQR ~ 0.0 

0.001 

IFUEL 

integer 

1 

1 -> no fuel is stored in the wing 

2 -> fuel is stored in the wing 


2 

CWMAN 

float 

1 

Design maneuver load factor 

6.25E+05 

1 

CF 



Shanley's const, for frame bending 


6.25E-05 














































































Fuselage 


Structural Concept 


Variable 

Type 

CKF 

float 

EC 

float 

KGC 

float 

KGW 

float 


imensio 



Units/Comment 


Description Unit! 


Frame stiffness coefficient. 


Power in approximation equation for 
buckling stability. 


Buckling coefficient for component general 
buckling of stiffener web panel. 


Buckling coefficient for component local 
buckling of web panel. 


KCON(T/B)| STRUCTURAL GEOMETRY CONCEPT 


2 I Simply stiffened shell, frames, sized for minimum weight in bucklin 


Z-stiffened shell, frames, best bucklin 


Z-stiffened shell, frames, buckling-minimum gage compromise 


Z-stiffened shell, frames, buckling-pressure compromise 


Truss-core sandwich, frames, best bucklin 


Truss-core sandwich, no frames, best bucklin 


Truss-core sandwich, no frames, buckling-minimum gage-pressure compromise 





0.505 


Default (747) 



Material Properties 


Variable 

Type 

FTS(T/B) 

float 

FCS(T/B) 

float 

ES(T/B) 

float 

EF(T/B) 

float 

DS(T/B) 

float 

DF(T/B) 

float 

TMG(T/B) 

float 

KDE 

float 

KDF 

float 

Geometric Parameters 

Variable 

Type 

CLBR1 

float 

ICYL 

integer 


imensio 



Dimensio 


1 



Description 


Tensile Strength on (top/bottom) 


Compressive Strength 


Young's Modulus for the shells 


Young's Modulus for the frames 


Density of shell material on (t/b) 


Density of frame material 


Minimum Gage thickness 


Knock-down-factor for modulus 


Knock-down-factor for strength 


Description 


Fuselage break point as a fraction 
of total fuselage length 


1 -> modeled with a mid-body 
cylinder. Else uses two power-la\ 
bodies back-to-back 


Units/Comment 



Default (747) 


58500 


54000 


1.07E+07 


1.07E+07 


0.101 


0.101 


0.071 


1 


1 


Units/Comment 

























































































Loads 

Variable 

Type 

Dimensio 

Description 

Units/Comment 

IMiHi IHHifMilM 

AXAC 

float 

1 

Axial Acceleration. 

s’ s 

0 

CM AN 

float 

1 

Weight fraction at maneuver. 


1 

ILOAD 

integer 

1 

1 -> analyze maneuver only. 






2 -> analyze maneuver and landing only. 

3 -> analyze bump, landing, and maneuver. 

3 

IT'lirt IWtit 

float 

12 

Fuselage gage pressure on (top/bot) 


13.65 

WFBUMP 

float 

1 

Weight fraction at bump. 


0.001 

WFLAND 

float 

1 

Weight fraction at landing. 


0.9 

Landing Gear 

Variable 

Type 

Dimensio 

Description 

Units/Comment 

Default (747) 

VSINK 

float 

1 

Design sink velocity at landing. 

10 

STROKE 

float 

1 

Stroke of landing gear. 

2.21 

CLRG1 

float 

1 

Length fraction of nose landing gear. 

0 1131 

CLRG2 

float 

1 

Length fraction of main landing gear 





measured as a fraction of total fuselage 

0.466 




length. 



WFGR1 

float 

1 

Weight fraction of nose landing gear. 

0.0047 

WFGR2 

float 

1 

Weight fraction of main landing gear. 

0.0398 

IGEAR 

integer 

1 

1 -> main landing gear located on fuselage 

2 -> main landing gear located on wing 

2 

GFRL 

float 

1 

Ratio of force taken by nose landing gear to 

0.001 




force taken by main gear at landing. 


CLRGW1 


1 

Position of wing gear as a fraction 

If only 1 wing gear. 

0.064 

CLRGW2 


1 

structural semispan. 

set CLRGW2 = 0.0 

0.1844 


Variable 

Type 

IT AIL 

integer 

Weights 

|| Variable 

Type 

jWTFF 

float 

[CBUM 

float 

)CLAN 

float 


Dimensio 



Description 


1 -> control surfaces mounted on tail. 

2 -> control surfaces mounted on wing 


Units/Comment 


ibesmbzhm 



Dimensio 


1 


Description 


Weight fraction of fuel. 


Weight fraction at bum 


Weight fraction at landing. 


Units/Comment 


Default (747) 


0.262 























































































Variable 

Type 

ISCHRENK 

integer 

ICOMND 

integer 

WGNO 

float 

SLFMB 

float 

WMIS 

float 

WSUR 

float 

wcw 

float 

WCA 

float 

NWING 

integer 


Dimensio 



Description I Units/Comment 


1 -> use Schrenk load distribution on wing 
else -> use trapezoidal distribution. 


1 -> print gross shell dimensions envelope 

2 -> print detailed shell geometry. 


Non-optimal factor for wing (including the 
secondary structure.) 


Static load factor for bumps. 


Volume component of secondary structure. 


Surface area component of secondary 
structure. 


Factor in weight equation for non-optimal 
weights. 


Factor in weight equation multiplying 
surface areas for non-optimal weights. 


Number of wing segments for analysis. 

































TABLE 8. KEY CONFIGURATION PARAMETERS FOR NINE TRANSPORT AIRCRAFT. 


ACSYNT INPUT PARAMETERS 

1. Geometry module 
Namelist WING 

"Variable"! 720 727 737 747 880 dc8 mdl 1 md83 1-1011 

SWEEP 35 32 25 37.17 33.5 30.6 35 24.16 35 

KSWEEP 222222222 

AR 6.958 7.67 8.21 6.96 7 7.52 7.5 9.62 6.98 

TAPER 0.333 0.2646 0.2197 0.2646 0.2494 0.1974 0.255 0.156 0.3 

TCROOT 0.1551 0.154 0.126 0.1794 0.08 0.1256 0.167 0.138 0.13 

TCTIP 0.0902 0.09 0.112 0.078 0.121 0.105 0.093 0.12 0.09 

ZROOT -1 -1 -0.25 -0.1 -1 -1 -0.79 -1 -1 

AREA 2460 1587 1005 5469 2000 2927 3648 1270 3590 

DIHED 336773633 

XWING 0.2963 0.376 0.35 0.249 0,301 0,302 0.218 0,468 0,359 

Namelist HTAIL 



Namelist VTAIL 



Namelist Fl)S 
















































































Namelist WPOD (inboard) 



Namelist WPOD (outboard) 



Namelist FPOD 



Namelist FPOD (third engine) 



2. Trajectory module 
Namelist TRDATA 






































































































3. Weights module 


Namelist OPTS 


Variable 

720 

727 

737 

747 

880 

dc8 

md 1 1 

md83 

1-1011 

ESH 

202000 

18202 

160000 

12759 

100800 

8165 

713000 

44290 

185000 

15158.2 

335000 

27058 

602500 

40955 

140000 

10340 

409000 

34797 


i 


PDCYL INPUT PARAMETERS 


Wing 

Geometric Parameters 


Variable 

720 

727 

737 

747 

880 

dc8 

mdl 1 

md83 

1-1011 

1STAMA 

2 

2 

2 

2 

2 

2 

2 

2 

2 

CS1 

0.1 

0.2125 

0.0724 

0.088 

0.1 

0.0818 

0.168 

0.181 

mm 

CS2 

0.27 

0.25 

0.238 

0.277 

0.264 

0.136 

0.2835 

0.271 

mESM 

Structural Concept 

Variable 

720 

727 

737 

747 

880 

dc8 

mdl 1 

md83 

1-1011 

CLAQR 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

IFUEL 

2 

2 

2 

2 

2 

2 

2 

2 

2 

CWMAN 

1 

1 

1 

1 

1 

1 

1 

1 

1 

CF 

6.25E-05 

6.25E-05 

6.25E-05 

6.25E-05 

6.25E-05 

6.25E-05 

6.25E-05 

6.25E-05 

6.25E-05 

Material 

Properties 

Variable 

720 

727 

737 

747 

880 

dc8 

mdl 1 

md83 

1-1011 

PS 

1 

1 

1 

1 

1 


1 

1 

1 

TMGW 

0.2 

0.2 

0.2 

0.2 

0.2 


0.2 

0.2 

0.2 

EFFW 

0.656 

0.656 

0.656 

0.656 

0.656 


0.656 

0.656 

0.656 

EFFC 

1.03 

1.03 

1.03 

1.03 

1.03 

1.03 

1.03 

1.03 

1.03 

ESW 

1.08E+07 

1.08E+07 

1.08E+07 

1.07E+07 

1.07E+07 

1.08E+07 

1.07E+07 

1.07E+07 

1.06E+07 

FCSW 

63500 

56000 

56000 

54000 

54000 

56000 

56000 

56000 

67000 

DSW 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

KDEW 

1 

1 

1 

1 

1 

1 

1 

1 

1 

KDFW 

1 

1 

1 

1 

1 

1 

1 

1 

1 


Fuselage 

Geometric Parameters 


Variable 

720 

727 

737 

747 

880 

dc8 

mdl 1 

md83 

1-1011 

CLBR1 

1.1 

1.1 

1.1 

1.1 

1.1 

1.1 

1.1 

1.1 

1.1 

ICYL 

1 

1 

1 

1 

1 

1 

1 

1 

1 



























































































Structural Concept 


Variable 



720 

727 

737 

747 

880 

dc8 

mdl 1 

md83 

1-1011 

5.24 

5.24 

5.24 

5.24 

5.24 

5.24 

5.24 

5.24 

5.24 

2.36 

2.36 

2.36 

2.36 

2.36 

2.36 

2.36 

2.36 

2.36 

0.368 

0.368 

0.368 

0.368 

0.368 

0.368 

0.368 

0.368 

0.368 

0.505 

0.505 

0.505 

0.505 

0.505 

0.505 

0.505 

0.505 

0.505 



Material Properties 

Variable 

720 

727 

737 

747 

880 

dc8 

mdl 1 

md83 

1-1011 

FTS(T/B) 

58500 

KgjiX 

58500 

58500 

58500 

64000 


|$i| 

58500 

FCS(T/B) 

54000 


54000 

54000 

54000 

39000 

54000 

ES(T/B) 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

EF(T/B) 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

DS(T/B) 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

DF(T/B) 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

TMG(T/B 

0.04 

0.04 

0.036 

0.071 

0.067 

0.05 

0.055 

0.055 

0.075 

KDE 

1 

1 

1 

1 

1 

1 

1 

1 

1 

KDF 

1 

1 

1 

1 

1 

1 

1 

1 

1 




Variable 


AXAC 

CMAN 

ILOAD 

PG(T/B) 

WFBUM 

WFLAN 



737 

747 

0 

0 

1 

1 

3 

3 

11.25 

13.65 

0.001 

0.001 

0.9 

0.9 



dc8 

mdl 1 

md83 

1-1011 

0 

0 

0 

0 

1 

1 

1 

1 

3 

3 

3 

3 

13.155 

11.5 

12.5 

12.6 

0.001 

0.001 

0.001 

0.001 

0.9 

0.9 

0.9 

0.9 



Variable 

720 

727 

737 

747 

880 

dc8 

mdl 1 

md83 

1-1011 

VSINK 

10 

10 

10 

10 

10 

10 

10 

10 

10 

STROKE 

1.67 

1.167 

1.167 

2.21 

1.33 

1.375 

1.9 

1.67 

2.17 

CLRG1 

0.133 

0.1306 

0.145 

0.1131 

0.095 

0.108 

0.141 


0.161 

CLRG2 

0.51 


0.5254 

0.466 

0.437 

0.499 

0.57 

1 

0.56 

WFGRl 

0.00389 

0.00725 

0.0052 

0.0047 

0.0057 

0.0311 

0.0031 

0.004 

0.005 

WFGR2 

0.036 

0.0738 

0.0382 

0.0398 

0.0322 

0.0742 

0.0058 

0.035 

0.044 

1GEAR 

2 

2 

2 

2 

2 

2 

2 

2 

2 

GFRL 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

CLRGW1 

0.1675 

0.1736 

0.1846 

0.064 

0.1574 

0.14 

0.2 

0.148 

0.232 

CLRGW2 

0 

0 

0 

0.1844 

0 

0 

0 

0 

0 




































































Factors 


















































PDCYLIN 


PS=1 . , 
EFFC=1 . 03 , 
KDEW= 1.0, 

ISTAMA=2 , 

CLAQR= .001, 

CKF=5 .24, 


TMGW= .02, 
ESW=10 . 7E06 , 
KDFW= 1.0, 

CS1=0 . 088 , 

IFUEL=2 , 

EC=2 .36, 


EFFW= . 656, 
FCSW=54000 . , 


CS2=0 .277, 
CWMAN=1 .0, 


KGC= .368, 


DSW=0 . 101, 


CF=6 . 25E-05 , 
KGW= .505, 


FTST = 4*58500 . , 8*0 . , 
FCST = 4*54000 . ,8*0 . , 
EST = 4*10. 70E06, 8*0. , 
EFT = 4*10 . 70E06, 8*0. , 
DST = 4* . 101, 8*0 . , 

DFT = 4* . 101, 8*0 . , 

TMGT = 4*. 07 1,8*0., 

KDE = 0.9, 


FTSB = 4*58500 . , 8*0 . , 
FCSB = 4*54000 . , 8*0 . , 
ESB = 4*10 . 70E06, 8*0 . , 
EFB = 4*10 .70E06, 8*0 . , 
DSB = 4* . 101, 8*0 . , 

DFB = 4* . 101, 8*0 . , 
TMGB = 4*. 07 1,8*0., 

KDF = 0.9, 


CLBR1=1 . 1, 

KCONT = 12*4, 

AXAC=0 . , 
CMAN=1 . 0 , 
WFBUMP=0 . 001, 

WTFF=0 .262, 

VSINK=10 . 0 , 
WFGRl = 0 . 0047 , 
CLRGWl = 0 . 064 , 

ITAIL=1, 

ISCHRENK=1 , 
WMIS=0 . , 
NWING=40 , 


ICYL = 1, 

KCONB = 12*4, 

CBUM=1 .0, 
ILOAD=3 , 
WFLAND=0 . 9 , 


STROKE=2 .21, 
WFGR2 = 0 . 0398 , 
CLRGW2 =0.1844, 


ICOMND=l, 
WSUR=0 . , 


CLAN=0 .791, 

PGB = 12*13.65, PGT = 12*13.65, 


CLRG1= . 1131 , 
IGEAR=2 , 


WGNO=l . 00 , 
WCW=1.0, 


CLRG2=0 .466, 
GFRL=0 .001, 


SLFMB=1 .2, 
WCA=0 .0, 


BEND 


Figure 16. PDCYLIN Namelist for 747-21P. 




^ DEG. DEG\ DEG. FT 3 LBS LB/FT3 

195.1006 114.493 41.4649 11.0904 0.247 0.179 0.078 42.98067 30.30835 40.200 28673.285 186806.00 8.504 
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Figure 17. PDCYL Wing Weight Output for 747-21 P. 







oooHrjmoohHinMn(N'fl(MH^HrMi£io\fMnO(MncDO'<j(na)inMrino^(Nco^oo'00(yia)(Dr'TT^^HntNoo\Hoo 

oooooo^<Nro'i3o^^(Nrnr^vDOOrH^inr^vc)fnir>cNja\ir>cD^<T\crvr^*--i*-Hvx)vx)CT>vx>tr»r^<Nr'U^rn^-iojr^'»s , inunir>Ln^rn(Noooo 

ooooooooooo«-icMm’^^chmmom® , a''^r^'H(N»-»*Hro*x>»-Hr^ ,, 3 , *HCDV£>^mc^t-H*Hoooooooooooooooooo 

ooooooooooooooooorHoi^vDiniDininvxivom^mrMrMiHrH^ooooooooooooooooooooooooo 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 




WAVGTO 




WAVGTO 




APPENDIX B: High-Altitude Study 


Description 

A study was made to estimate the wing weight of a scaled version of an existing 
propeller-driven high-altitude drone aircraft. This aircraft, termed the Strato7, is modeled 
as an enlarged version of the existing Perseus-a3. PDCYL was used to validate the wing 
weight estimation returned by ACSYNT. 

The wing of the Strato7 incorporates a single hollow, cylindrical carbon- 
fiber/epoxy spar placed at the leading edge. The strength of the cover is assumed 
negligible. No fuel is carried in the wing, while propulsion and landing gear are mounted 
on the fuselage. The layout of the Strato7 is shown in Figure 20. 

Input 


Fuselage weight estimation is not considered for the Strato7. An example of the 
ACSYNT input for the Strato7 wing weight estimation is shown in Figure 21 . The 
corresponding PDCYLIN namelist for the case where the ratio of structural chord to total 
chord is 0.2 is shown in Figure 22. 

Output 

Wing weight as a function of the ratio of structural chord to total chord is shown 
in Figure 23. The wing weight estimated by ACSYNT is 789 pounds. PDCYL matches 
this wing weight when the ratio of structural chord to total chord is approximately 0.25. 
Non-optimum weight was not considered in this analysis. In order to estimate non- 
optimum weight, non-optimum factors would need to be recomputed for this type of 
aircraft. 



Figure 20. Strato7 Configuration. 
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Figure 21. ACSYNT Input for Strato7. 



0 , 

0 , 


SEND 


o o 


EPDCYLIN 


PS=1 . , 

TMGW= .05, 

EFFC=1 .108, 

ESW=12 . 9E06 , 

KDEW=1 .0, 

KDFW=1 .0, 

ISTAMA=2 , 

CS1=0.01, 

CLAQR= . 001 , 

IFUEL=1 , 

CKF=5 . 24 , 

o 

o 

CN 

II 

U 

W 

FTST = 4*58500. 

, 8*0. , 

FCST = 4*54000. 

,8*0. , 

EST = 4*10 .70E06, 8*0 . , 

EFT = 4*30. 0E06, 8*0. , 

DST = 4*. 101, 8* 

0. , 

DFT = 4* .292, 8* 

0. , 

TMGT = 4* .03, 8* 

0. , 

KDE=0 . 9, 


CLBRl=l . 1 , 

ICYL = 1, 

KCONT = 12*4, 

KCONB = 12*4 


EFFW= . 605 , 

FCSW=75000 . , DSW=0 .058, 

CS2=0 .75, 

CWMAN=1 . 0 , CF=6 . 2 5E-05 , 

KGC= . 368 , KGW= . 505 , 

FTSB = 4*58500. , 8*0. , 

FCSB = 4*54000 . , 8*0 . , 

ESB = 4*10 . 70E06 , 8*0 . , 

EFB = 4*30 . 0E06 , 8*0 . , 

DSB = 4*. 101, 8*0. , 

DFB = 4* .292, 8*0 . , 

TMGB = 4* . 03, 8*0 . , 

KDF=0. 8, 


AXAC=0 . . 
CMAN=1 . 0 , 
WFBUMP=0 . 001, 


CBUM=1 .0, 
ILOAD=3 , 
WFLAND=0 .9, 


CLAN=0 . 93 , 

PGB = 12*11.5, PGT = 12*11.5, 


WTFF=0 . 07, 


VSINK=10 .0, 
WFGR1=0 . 0031 , 
CLRGW1=0 . 20 , 

ITAIL=1, 


STR0KE=1 . 0 , 
WFGR2=0 .0058, 
CLRGW2 = 0.0, 


CLRG1= . 395 , 
IGEAR=1 , 


CLRG2=0 .5, 
GFRL=0 . 001, 


ISCHRENK=1 , 
WMIS=0 . , 
NWING=40, 


IC0MND=1 , 
WSUR=0 . , 


WGN0=1 . 00 , 
WCW=1 .0, 


SLFMB=1 .2, 
WCA= 0.0, 


$END 


Figure 22. PDCYLIN Namelist Input for Strato7. 



Figure 23. Strato7 Wing Weight as a Function of Structural Box Length. 


